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Abstract 

Using the Bethe ansatz we obtain in a determinant form the exact 
solution of the master equation for the conditional probabilities of the 
totally asymmetric exclusion process with particle-dependent hopping 
rates on Z. From this we derive a determinant expression for the time- 
integrated current for a step-function initial state. 

1 Introduction 

In the usual one-dimensional totally asymmetric simple exclusion process 
(TASEP) particles hop randomly in continous time on the integer lattice Z. 
A hopping event occurs independently for each particle after an exponen- 
tially distributed random time with parameter 1 to its right neighboringsite, 
provided this site is empty; otherwise the attempted move is rejected jH 1^ ■ 
This is the exclusion interaction between biased random walkers, which phys- 
ically models a short-range hard-core repulsion of identical driven diffusive 
particles. Here we consider the case of distinct particles where each particle 
i has its own intrinsic hopping rate fj. 
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This model has been studied earher by other authors. To our knowledge it 
was first introduced (with a random choice of rates w drawn from some given 
distribution) by Benjamini et al. a who proved the existence of a product 
invariant distribution (as seen from a tagged particle) above some critical 
density p* of particles. They also proved for this case the hydrodynamic limit 
of an associated zero-range process (ZRP) which is obtained from the TASEP 
by identifying particles with the sites of a new 1-D lattice and the interparticle 
distance (number of empty sites between particles i,z + 1) as occupation 
number at site i of that lattice. Subsequently it was shown for periodic 
boundary conditions that below the critical density p* a macroscopic traffic 
jam, trailing the slowest particle, develops The macroscopic dynamical 

signature of the transition under Eulerian scaling is a discontinuity in the 
coarse-grained density at the front of a rarefaction wave that develops out 
of an initial step function profile The stationary distribution of an open 
system where particles are injected and extracted was obtained numerically 
|7|. Somewhat surprisingly, the theory of boundary- induced phase transitions 
developed by Kolomeisky et al. [1] accounts for the stationary phase diagram 
of the system as a function of the effective boundary densities. 

The jamming transition occurring in this model is a classical analog of 
Bose-Einstein condensation and one of the motivations to study this system 
in more detail. An analogous transition occurs also if the rates are identical 
for each particle, but dependent on the lattice distance to the next particle, 
see joj for a recent review (in terms of the ZRP) and 1^, 11, ^ 3 JJ,E| 



for the current developments. A second motivation comes from the recent 
interest in particle systems with several conservation laws jl^. By choosing 
the Vi in some way out of p distinct values the model indeed becomes a particle 
system with p conservation laws where, unlike in a similar class of models 
introduced by Karimipour 17, 1^ 3' overtaking of particles is permitted. 
Exact results for such multi-component systems are scarce. Our contribution 
aims at providing some tools for future exact and rigorous analysis. We 
first solve the master equation for an A^-particle system defined on Z with 
the Bethe ansatz generalizing the approach taken by one of us earlier j2ol |. 
This yields a determinant formula for the conditional probabilities which 
we then use to derive an exact expression for the distribution of the time- 
integrated current for the step initial configuration envisaged in Our 
resu lt g eneralizes that obtained recently for the TASEP with constant rates 



21, 2 
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2 Solution of the master equation 



We consider initial states such that there is no particle beyond a lattice 
point Xmax < oo at time t = 0. Since the process is totally asymmetric and 
defined on Z the rightmost particles are not infiuenced by particles to their 
left. Without loss of generality we may therefore consider A^-particle initial 
states where the particles are labelled from 1 through N from left to right. 
Therefore at any finite time time t > a configuration of the system can be 
represented by an increasing sequence of integer numbers {xi,X2, - ■■ ,xn), 
Xi denoting the position of the ith particle with hopping rate f j. We assume 
< f j < fmax where fmax < oo sets the microscopic time scale. We remark 
that allowing for fmax = cxd is equivalent to considering particles covering 



more than 1 (neighboring) site |23l . |2J, |25|, |26|, |27|, |28[ . Such a system can 



be mapped to an exclusion process with standard particles (of lattice size 
1) by an appropriate coordinate shift. The case f j = is trivial in so far 
as this amounts to cutting the lattice Z into finite segments in which after 
some time all particles pile up at the right boundary of the segment and the 
dynamics freeze. 

In case of equal hopping rates the exchange of particles is irrelevant as 
they are indistinguishable. However, in this generalized version one has to 
define the corresponding exchange rules separately since the particles have 
their own identity. In our model particles cannot exchange so one can in- 
terpret it as a one-lane traffic model with cars having their own preferred 
speed. 

The master equation for the probability of finding particles at sites Xi is 
the following (the t-dependence is dropped for simplicity of notation): 

^P(X1, X2, • ■ ■ , Xn) = ViP{xi - 1, X2, ■ ■ ■ Xn) 

+ f2(l - 5^2_l,a;i)-P(Xl, X2-1,..., Xn) 

h Vn{1 - 5^i^-l,XN-i)PiXl^ X2 - 1, . . . , Xjv - 1) 

- {Vi{l - 42-l,xi) + V2{1 - 5z3-l,a;2) H \- Vn)P{Xi,X2, ■ ■ .,Xn), (1) 

where P{xi, X2:r-- , Xn) is the probability of the configuration (xi, X2, ■ ■ ■ , Xn) 
In the spirit of j20| we extend the range of definition of P(xi, X2, - ■ ■ ,X]\i) from 
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the physical domain < Xj+i to Z . By requiring the boundary condition 

f fc+lP(xi, ■ ■ ■ , Xk-l, Xfc, Xk, Xk+2, ■ ■ ■ , Xn) 

= VkP{xi, ■ ■ ■ , Xk-l, Xk, Xk + 1, Xk+2, ■ ■ ■ , Xn) (2) 
to be vahd for alH > the master equation reduces to the simple form 

N 

, Xn) = ^ ] ViP^Xi, . . . , Xi-i, Xi — 1, ■ ■ • , Xn) 

i=l 
N 

- ^ ViP{xi, X2,..., Xn) = -HP{xi, X2, ■ ■ ■ , Xn)- (3) 

i=l 

The second equality defines the linear operator H generating the time evolu- 
tion of the probability distribution. For given initial configuration [yi, y2, - ■ ■ , Vn) 
this quantity is the conditional probability P{xi, ■ ■ ■ , xnIvi, ■ ■ ■ , Vn] t)- 

It is obvious that (jHl) can be solved by Fourier transformation with a "mo- 
mentum" variable ki associated with each particle coordinate Xi. However, a 
straightforward Fourier ansatz does not satisfy the boundary condition Q. 
This problem can be resolved by noting that there is some freedom by choos- 
ing a linear combination of Fourier transforms over all permutations of the 
dummy momentum variables ki with complex coefficients (wave amplitudes) 
that depend on the ki. It is the essence of the Bethe ansatz to take this 
approach and to implement it by factorizing the wave amplitudes for the 
A^-particle problem into amplitudes for the two-particle problem. It is done 
such that for each permutation of the momentum variables the iV-particle 
amplitude picks up the corresponding two-particle amplitude as extra factor. 

In the present case this strategy also turns out to be successful, but re- 
quires some modification due to the occurrence of the different hopping rates 
f j. In order to obtain the eigenvalues and eigenfunctions of the "Hamilto- 
nian" H with the boundary condition Q we write 

^Pk{xi, ■ ■ ■ , Xn) = -E]^Pk{xi, ■ ■ ■ , Xn)- (4) 

where k represents the conjugate momenta ki, . . . , kN- Modifying the Bethe 
ansatz to account for the Vi leads to eigenfunctions of the form [l^ 

TV / ^ \ 

Pfci,-,fcjv(^i) ■ ■ • )^Af) = JJ ^^i'^' 2Z -S"^ exp ii^Xjk^^j) j , (5) 

i=l TT \ j = l / 



— P{Xi,X2, 
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where the sum is taken over all permutations vr of (1, 2, 3, ■ ■ ■ , A^). The above 
eigenfunctions correspond to the eigenvalues 



N 



fell'" jfciV 



(6) 



With the inverse Fourier transformation 



N 

n 

i=l 



dk- 

Z7l 



(7) 



it is trivial to show that this ansatz satisfies (0) with arbitrary wave ampli- 
tudes S'^(k) . 

On the other hand from the boundary condition Q it follows that 



^ 1 _ v^e^'-ik+i) 



(8) 



where the operator changes the kth and k+ 1st element of vr and S(^i^2,...,n) 
is an arbitrary function defined by the initial condition to be satisfied at 
t = 0, see next section. 

In order to obtain the A^-particle amplitude 5*,^ for a given permutation 
71 one has to decompose this permutation into a sequence of nearest neigh- 
bour exchanges and take the product of factors appearing in (jH)). Although 
this decomposition is not unique the resulting S',r is well-defined if certain 
relations, known as Yang-Baxter equations, are satisfied. Instead of checking 
this explicitly one can easily show that the r.h.s. of (0) with an S'tt satisfying 
(jHI) can be written in the following determinant form 



Pkl,-,kN{Xl; 

S{i,...,N)Ylvi^^ det 



«jvv-i) ,Xn) 
N 



i=l 



n— 1 m—l 



■ (9) 



m,,n 



This verifies uniqueness by construction and therefore also implies the Yang- 
Baxter equations. 
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3 Conditional probability 



In order to obtain the conditional probability for the initial configuration 
(yi, ■ ■ ■ ,yN) we require 

P(x; t = 0) = 5x,y. (10) 
In the following we show that the choice 

(N \ ^ 

-ij^y^'^AH^i^' (11) 
j=i J i=i 

with an appropriate (see later) definition of the contours of integration in the 
complex ki planes (|7j) indeed reproduces the required initial condition. 

In order to write the conditional probability in a compact form we define 
the functions 

^ „ l-l k-l 

1=1 i=l 

Here, as throughout this paper, the empty product is defined as unity. The 
integral has to be taken along a circle of radius e around the origin of the 
complex plane. Notice that if all the rates are less than a given number fmax 
(as we assume here) then the integral can be taken just as well along the 
circle with radius 1/fmax — 0. Now we exploit the determinant expression ^ 
to obtain 

Theorem 1. (Conditional probability) LetP(x|y;t) be the conditional prob- 
ability for the TASEP with N particles on the infinite one- dimensional lattice 
with right hopping rates fi,f2, ■ ■ ■ iVn (from left to right), starting at sites 
yi,y2,--- ,yN at time t = 0. The conditional probability of finding these 
particles on sites Xi,X2, ■ ■ ■ ,xn at time t can be written as the determinant 

N 

P(x|y;t) = J](e-*^>i;,^-^0 

i=l 



i^i,i(a;i 
i^2,i(a;i 


-y2,t) 


Fl,2{x2 
F2,2ix2 


-yi,^) ■ 
-y2,t) ■ 


Fi^n^xn 

F2,n{xn 


-yi^t) 
-y2,i) 


Fn,i{xi 


-yN,t) 


FN,2ix2 


-yN,t) ■ 


■ Fj^^]\r{x]\f 


- yN,t) 



(13) 
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Remark. The expression (jl3p is a generalization of the result found in [2(|] for 
the case of equal hopping rates Vi = 1. The functions Fk^i{x,t) then reduce 
to e^Fi_k{,x,t) of that work. 



Proof. One has to show that P(x|y;t) satisfies (j21), © and (jTU)). Using the 
specific form of Pk given by Q and (fTT|) the inverse Fourier transformation 
© leads to 



N 

n 



2^' 



e ■^'=*det 



n-l 

.i=i 



m— 1 



(14) 

This expression satisfies ((21) and (jHI) by construction, independently of the 
contour of integration. Notice that the different rows of the determinant 
contain different k variables which enables us to perform the integration 
over ki separately on each element of row i. After changing the variable k to 
z = e^^ (|T^ indeed reduces to the r.h.s. of (fT!?|) . By adopting the argument of 
j2fll | one can easily show that the specific choice of the contour of integration 
in (O indeed yields □ 

In the following F{x) has to be understood as F{x,t). 
Corollary. Consider the probability 

P(z|y;t) = Prob(xi > Zi,X2 > ^2, ■ ■ ■ ,Xn > ZN\y;t) 

that all N particles, starting on y, have reached "at least" the configuration 
z where zi < Z2 < ■ ■ ■ < . This quantity can be written as a determinant 



N 



P(z|y;t) = J](e-*^'t;/-^') 



i=l 



X 



Fi,2{zi-yi) Pl,3(^2-l/l) 

^2,2(2:1 -^2) ^^2,3(^2 - 1/2) 
Fn,2{zi - yn) Fn,3{^2 - Vn) 



Fi,n+i{zn - yi) 
F2,n+i{zn - 1/2) 

Fn,n+i{zn — yN) 



(15) 
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Proof. One has to write -P(z|y; t) in the form 

OO — 1 X4 — 1X3 — 1X2 — 1 

p(z|y;t)= J2 E ■■■EE E^(^'y;^) (i^) 

Xn=ZN Xm-1=ZN-1 X3=23 X2=22 Xi=2i 

then use (fT^. Applying the identities 

X2 

^ = - t;;"^+^Ffc,;+i(a;2 + 1) (17) 

x=x\ 

X2 

E = Vk-i'^'Fu-i^iix,) - Vk-i^'+'Fk-i,i{x2 + 1) (18) 



X = Xl 



and standard determinant-manipulations then leads to (fT3j) . The identities 
()17ll8p follow by straightforward computation from the definitions ()12|). □ 

4 Current distribution 

We consider now the step function initial state where the lattice is fully 
occupied up to the lattice site x = and empty for x > 0. The time 
evolution of the mean density (under Eulerian scaling) has been obtained 
in for a large class of random distributions of rates fj. An interesting 
problem are the fluctuations of the current which were first obtained in i29[ 
for the homogeneous system with Vi = 1. As a step towards characterizing 
the long-time properties of these fluctuations (which we do not consider) we 
derive here for flnite times an exact expression for the probability that the 
time-integrated current across a given bond (i.e. the number of particles that 
have crossed that bond up to time t) equals some number N. This quantity 
is related to the distribution of the distance travelled by the flrst particle 
with the initial condition where N particles sit on consecutive sites. More 
precisely, 

Theorem 2. (Current distribution) choosing as initial sites the set yN = 
(1 — A^, ■ ■ ■ ,0) the probability that all particles have crossed the bond x — l,x 
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up to time t can be written as 



N 



QNix,t) = l[{e-'^^v,^-'''~') 



i=l 



Fi,2(x + Ar-1) Fi,3(x + Ar) 

F2,2{x + N-2) F2,3{X + N-1) 



Fi,;v+i(x + 2iV-2) 
F2,N+i{x + 2N - 3) 



(19) 



X 




F^^N+^{x + N - 1) 



Remark. In the case of the initial condition yoo when the lattice is fully 
occupied on the left and empty on the right of the origin Qjq{x,t) is the 
probability that at least N particles have crossed the bond a; — 1, x up to 
time t. On the other hand this is equal to the probability that the A^th 
particle (from the right) reached at least the site x up to time t. 

Proof. Notice that Qjq{x,t) = P{x,x + 1, ■ ■ ■ ,x + N — l|yN;t). By using 
(fT^ one arrives at (fT^ . □ 

Remark. The determinant appearing in ()19|1 has several equivalent forms. 
One can derive these by using the identities 



which are special cases of ()17|18p corresponding to Xi = X2. An example is 
the following 



Fk,i{x) = Fk,i+i{x) - viFk,i+i{x + 1) 
Fk,i{x) = Fk-i,i{x) - Vk-iFk-i^i{x + 1), 



(20) 
(21) 



TV 



QN{x,t) = \{{e-''''^vr''-') 



i=l 



Fi,2(a; + iV-l) Fi,3(a; + Ar) 
Fi^2{x + N-2) Fi,3(x + Ar-1) 



Fi,^+i(x + 2iV-2) 
Fi,^+i(x + 2iV-3) 



(22) 



X 




Fi,z{x + l) 



i^i,TV+i(x + iV-l) 
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Another equivalent expression is 



N 



QN{x,t) = l[{e-'^^vr''~') 



i=l 



X 



F,,N+i{x + N - 2) F^^N+iix + N - 1) ■■■ Fi,^+i(x + 2Ar-3) 
Fi,N+i{x) Fi^N+i{x + l) ■■■ F,^r^+i{x + N - 1) 



(23) 



5 Final remarks 



The study of current fluctuations provides insight into universal properties 
of far-from-equilibrium systems and therefore has attracted a considerable 
amount of attention in recent years Q Q H H H H 111 H H 



l36l l37| . There is, however, to our knowledge no work on current fluctua- 
tions in particle systems with more than one conservation law. The result 
(|23|) is a first step in this direction. The form of (|23p indicates that the 
distribution QN{x,t) is independent of the order of the hopping rates of 
particles (e.g. vi = a,V2 = b,V3 = c would lead to the same QN{x,t) as 
f 1 = c, f 2 = ct, ^^3 = b) since these hopping rates enter symmetrically in 
Fi^N+i{x,t). This at first sight perhaps counterintuitive property becomes 
obvious in the growth presentation of the exclusion process P| and also in the 
last passage percolation picture 0, 11^ where the hopping rates Vi represent 
line defects with a different distribution of random energies. For the statisti- 
cal weight of paths crossing these defects it is clearly irrelevant in which order 
they occur. In the interpretation of the TASEP with particle-dependent rates 
as an A^-component lattice gas this means that the nonergodicity of the pro- 
cess implied in the absence of passing is irrelevant for the fluctuations of the 
total current. 

A next step would be to study the long time behaviour of the current 
fluctuations. This was calculated recently for the homogeneous ASEP 21, 2^ 
HHElllJ and the polynuclear growth model (PNG) 0,1310. It was 
found that the current fluctuations in the long time limit are characterized by 
universal scaling forms, which are given explicitly in terms of Tracy- Widom 
distributions known from random matrix theory. It would be of interest to 
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see whether this universahty is preserved in the case of particle- wise disorder. 

Recently Priezzhev 3q] has extended the Bethe ansatz solution for the 
usual TASEP on Z to periodic boundary conditions. It would be interesting 
to generalize this analysis to the present case. It would also be interesting 
to investigate whether determinant solutions can be obtained for integrable 
models with passing [l7, 3, 3 2^- 
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